Dynamics of Localization Phenomena for Hardcore Bosons in Optical Lattices 
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We investigate the behavior of ultracold bosons in optical lattices with a disorder potential gener- 
ated via a secondary species frozen in random configurations. The statistics of disorder is associated 
with the physical state in which the secondary species is prepared. The resulting random potential, 
albeit displaying algebraic correlations, is found to lead to localization of all single-particle states. 
We then investigate the real-time dynamics of localization for a hardcore gas of mobile bosons which 
are brought into sudden interaction with the random potential. Regardless of their initial state and 
for any disorder strength, the mobile particles are found to reach a steady state characterized by 
exponentially decaying off-diagonal correlations and by the absence of quasi-condensation; when the 
mobile particles are initially confined in a tight trap and then released in the disorder potential, their 
expansion is stopped and the steady state is exponentially localized in real space, clearly revealing 
Anderson localization. 

PACS numbers: 03.75.Lm, 03.75.Mn, 64.60.My, 72.15.Rn 



I. INTRODUCTION 

The adiabatic loading of ultracold atoms and molecules 
in optical lattices represents a formidable opportunity 
to engineer strongly correlated states of quantum many- 
body systems with unprecedented contro l 1 ' 2 ' 3 ' . On the 
one hand, known fundamental models for the physics of 
solid state systems can be literally implemented in op- 
tical lattices^, and the experimental detection of their 
ground state or thermal equilibrium properties has the 
potential of responding to the lack of theoretical results 
due to fundamental limitations in the calculations, as 
those generally reported in fermionic systems or in frus- 
trated quantum magnets. On the other hand, optical 
lattice systems offer the further advantage of control- 
ling the Hamiltonian parameters in real time, and this 
enriches the range of correlated phases that can be im- 
plemented, if one can guide the evolution of the system 
towards an off-equilibrium state which is not necessarily 
an eigenstate of a known Hamiltonian. 

An intense activity has recently been fo- 
cused on the experimental implementation of 
disorder potentials in systems of trapped ultra- 
cold bosons, both experimentally- 7 -' 8 -' ^ 10 ' 11 ' 12 and 
theoreticallyi 3 ^^ 1 ^^^^^. The introduc- 
tion of tunable randomness in the system offers the 
possibility of realizing Anderson localization of coherent 
matter waves of weakly interacting boson o 8 ' 9 ' 10 and 
the opportunity of investigating the interplay between 
strong localization and strong interaction, e.g. in a deep 
optical latticeii. The fundamental model describing 
this rich phenomenology in presence of a lattice is the 
Bose-Hubbard model in a random potential 24 , where, 
beside the conventional Mott insulating and superfluid 
phases, a Bose-glass phase appears, either associated 
with the fragmentation of weakly repulsive bosons into 
exponentially localized states, or with the localization of 
collective gapless modes for strong repulsion. 

The disorder potential has so far been realized opti- 



cally, either through laser speckles 7 -' 8 -' 9 -' 1 ^ with or without 
an optical lattice, or, in optical lattices, through a sec- 
ondary incommensurate standing wavei 1 -. In the case of 
speckles, the typical length scale associated with the dis- 
order potential is quite extended with respect to the cor- 
relation length of the bosons, so that classical trapping 
rather than quantum localization is responsible for the 
observed suppression of transport propertie o 12 ' 25 . On the 
other hand, the incommensurate superlattice of Ref. [ll| 
realizes a potential which is strongly fluctuating over the 
distance of a few lattice sites, but its psewdo-disordered 
nature requires to take also into account gapped insulat- 
ing phases at incommensurate fillings, competing with 
the Bose glass 2 ^. 

An alternative proposal to create disorder in optical 
lattices involves the repulsive interaction with a sec- 
ondary species of particle a 27 ' 28 . If the two species cor- 
respond to two different hyperfine states of the atoms, 
the use of state-dependent optical lattice a 29 ' 30 ' 31 allows 
one to first decouple the two species, and then to selec- 
tively suppress the hopping of one of them, freezing it in 
a given quantum state |$ f ). 

If the two species are subsequently brought into inter- 
action, the mobile one evolves in a quantum superposi- 
tion of all possible realizations of the random potential 
associated with the Fock components of |$ f ), and expec- 
tation values over the evolved state are therefore auto- 
matically averaged over the disorder distribution^ 8 -. In 
this proposal not only is the disorder potential strongly 
fluctuating over the length scale of a few lattice sites, but 
in principle one can also vary its statistics by preparing 
the frozen bosons in different states |4> ). 

In this paper we investigate the above proposal in de- 
tail in the exactly solvable case of a one-dimensional 
gas of hardcore bosons on a lattic o 32 ' 33 . We consider 
both species of bosons (the mobile one and the frozen 
one) to be hardcore repulsive, which is not only ex- 
perimentally feasibl o 33 ' 34 , but it has three fundamental 
theoretical advantages: 1) the state |$ f ) of the frozen 
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bosons can be obtained exactly from the Hamiltonian 
of the system before freezing, and disorder averaging 
can therefore be accurately performed; 2) Jordan- Wigner 
diagonalization 32 allows to calculate the exact real-time 
evolution of the mobile particle a 35 i 36 ' 37 after they are 
brought into interaction with the frozen ones; 3) and, 
most importantly from a conceptual point of view, lo- 
calization phenomena of hardcore bosons are perfectly 
understood in terms of Anderson localization of non- 
interacting fermions. While the observation of localiza- 
tion for interacting bosons is limited by screening of the 
disorder potential and the reduction of the healing length 
due to the interactional 9 - ' 10 ' 11 ' 12 , many-body effects en- 
ter the system of non-interacting fermions only through 
Fermi statistics. 

In particular we focus on the case in which the frozen 
bosons are initially prepared in the superfluid ground 
state | <E> ) of the hardcore-boson Hamiltonian at half fill- 
ing with periodic boundary conditions. The resulting 
disorder potential has the structure of a bimodal random 
on-site energy, with algebraically decaying correlations, 
but also with a very rich Fourier spectrum dominated by 
short-wavelength components. Most importantly, such a 
potential is found to lead to Anderson localization of all 
single-particle eigenstates (apart from possibly a set of 
zero measure) as expected for uncorrelated disorde r 38 ' 39 . 
The evolution of the mobile bosons in such a potential is 
found to invariably lead to a disordered steady state with 
exponentially decaying correlations and absence of quasi- 
condensation for a wide variety of realistic initial condi- 
tions. In particular phenomena of quasi-condensation in 
finite-momentum states 3 ^ and fermionizatior* 3 -^, reported 
upon expansion of the hardcore bosons from a Mott- 
insulating state and a superfluid state respectively, are 
completely washed out by the disorder potential. Hence 
we can conclude that this setup allows for a robust im- 
plementation of a localized state as the off-equilibrium 
steady state of the system evolution. 

This paper is structured as follows: SectionHTldescribes 
the system of two bosonic species, the exact diagonal- 
ization method for the study of real-time evolution, the 
Monte Carlo sampling of the disorder distribution, and 
the main features of the random potential; Section [Ml is 
devoted to the study of localization of the single-particle 
eigenstates in the random potential; Section IIVI investi- 
gates the evolution of the mobile bosons after interaction 
with the frozen ones when both species are prepared on 
a ring; finally, Section [V] is dedicated to the study of the 
expansion of the mobile bosons in the random potential, 
starting from different initial confined states. 



II. SYSTEM AND METHOD 

In this section we present the system of two bosonic 
species which is used to study the effect of a disorder 
potential (subsection III A[) . We then briefly describe the 
numerical procedure to exactly treat the equilibrium and 



out-of-equilibrium properties of the system (subsection 
III Bp , and the sampling of the disorder distribution (sub- 
section HTC]). Finally, we examine the nature of the cor- 
relations in the disorder potential created by a frozen 
species of hardcore bosons in subsection III Dl 



A. Hamiltonian Dynamics 

The full system of two trapped interacting bosonic 
species in a one-dimensional optical lattice is described 
by the Hamiltonian 
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where Tlo is the Hamiltonian of the bosons which will 
remain mobile, 
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Hf is the Hamiltonian for the bosons to be frozen, 
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, (3) 



and Hint is the interaction Hamiltonian 



(4) 



Here L is the number of sites of the system, a] and 
cii are boson creation and annihilation operators and 
rii = a\ai is the number operator for site i. Symbols 
with the superscript f are the corresponding operators 
for the frozen bosons that create the disorder potential. 
In Eqs. we also consider the possibility of both 

species being confined by a parabolic trap with different 
trapping strengths V, V 1 . 

At times t < the two species of bosons are not inter- 
acting with each other (W = 0), and they are prepared 
in the factorized ground state | \P) = | $ f ) (g> | $) of their 
respective Hamiltonian with fixed numbers of particles 
N and N l . The ground state |<i> f ) of the frozen particles 
can be decomposed in the Fock basis 



l* f > = E c (K})IK}>. 
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where the sum extends over all Fock states, 



Making use of the Jordan- Wigner transformatio: 



,32 



{ n i} = n n n 2i 



(6) 



At some time t < the frozen bosons are made immo- 
bile (J f = 0), and subsequently at t = the interaction 
between the two species is turned on (W > 0). Further- 
more, releasing the mobile bosons from their trap (V = 0) 
allows us to study their expansion properties. 

The time evolution of the initially prepared state for 
t > is described by 



|vl/(<))^ e -^^ c ({n,f})|{n l f })®|$) 
{"0 
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(7) 



where the _ coefficients ctfnj}, t) = e^W* c ({nj}) 
have acquired a phase factor which will become irrele- 
vant, and the state 

|$ (t, {n[})) = exp {-i [Ho + H mt {{n{})] t} |$) (8) 

represents the time evolution of the initial state of the 
mobile bosons interacting with a single Fock state {n?}) 
of the frozen bosons, which determines the static external 
potential. We have used the property that, for J = 0, 

[H{, Ho + Hint] = 0. 

Hence, equation (J7J) describes the parallel time evo- 
lution of the mobile bosons in a quantum superposi- 
tion of different realizations of the disorder potential 
Vi = Vn\, each appearing with a probability P ({n^}) = 

l c ({ n i}) I 2 - Remarkably, the time evolution of the ex- 
pectation value of an operator A acting only on the mo- 
bile bosons is automatically averaged over the disorder 
statistics^: 

(*(*)|A|*(i)> = 

£ \c({nl})\ 2 ($(t,{nl})\A\$(t,{nl})). (9) 

{«0 



B. Hardcore limit and Jordan- Wigner 
Transformation 



From here onwards we will restrict ourselves to the ex- 
actly solvable case of hardcore bosons, obtained in the 
limit U, U l — > oo for filling smaller than or equal to 
one. It is convenient to incorporate the hardcore con- 
straint directly in the operator algebra, passing to hard- 
core boson operators which anfzcommute on the same 
site, {a,, al} = 1, {a^,a^} = 0, and = 1, 
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(10) 



one can map the hardcore bosons operators onto spinlcss 
fermion operators fj and obeying the same Hamil- 
tonian as the bosonic one, apart from a boundary term 
depending on the number of particles in the case of pe- 
riodic/antiperiodic boundary conditions. The fermionic 
problem is exactly solvable, and its eigenstates can be 
written in the following form 



N 
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(ii) 



where the N columns of the matrix P nm — {P} nm repre- 
sent the first N single-particle eigenstates. Following the 
recipe of Ref. l35l . from the matrix P nm one can efficiently 
calculate the one-particle density matrix (OPDM) 



Pi j = (4 a j)> 
and hence the momentum distribution 

L 

L 
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(12) 



(13) 
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which represents a fundamental observable in trapped 
atomic systems. At a more fundamental level, the knowl- 
edge of the eigenvalues A,, of the OPDM, associated to 
eigenvectors $ also known as natural orbitals (NO), al- 
lows one to rigorously study condensation phenomena 
through the scaling of the maximum eigenvalue Ao with 
the particle numbe r 35 ' 40 . In absence of an external po- 
tential the OPDM of the hardcore bosons decays al- 
gebraically as pij ~ \i — j\ a , where a = 0.5, signal- 
ing off-diagonal quasi- long-range order in the system 41 . 
Correspondingly the occupation of the k — momen- 
tum state, which coincides with the natural orbital with 
largest eigenvalue for a translationally invariant system, 
scales as nk=o ~ VN with the particle number TV for 
any fixed density n = N/L < 1, namely it exhibits quasi- 
condensation. Remarkably, quasi-long-range order and 
quasi-condensation (in the form of a y/N scaling of the 
largest eigenvalue Ao of the OPDM) survive also in pres- 
ence of a trapping potential V(i — io) when the particle 
number is increased and correspondigly the trap strength 
V is decreased so that the characteristic density in the 
trap 



"= v '7 



(14) 



0. 



is kept constant and smaller than a critical value p c (w 
2.6 for a — 2) to avoid formation of a Mott plateau in 
the trap center—. 



4 



Finally, the exact solution of the fermionic Hamilto- 
nian allows one to calculate the real-time evolution of 
the fermionic wavefunction, and in particular of the P 
matrix as 



P(t) 



iHti 



(15) 



with the single-particle time evolution operator given by 



(e-* Ht L = <0|/, 



-mt a 



f Jo). 



(16) 



Making use of this approach, Refs. l36ll37l have shown that 
quasi-condensation is a robust feature of the system af- 
ter expansion starting from an initially trapped quasi- 
condensed state, and it is even dynamically recovered 
when the initial state before expansion is a fully incoher- 
ent Mott insulator state. 



C. Disorder Averaging 

Eq. ((9|) shows that, ideally, the unitary evolution of 
the system explores all possible realizations of the disor- 
dered potential at once. Numerical calculations based on 
a matrix-product-state representation of the system state 
also enjoy this feature of "quantum parallelism" of the 
Hamiltonian evolution 28 by treating the disorder poten- 
tial as a quantum variable in the system. In this paper we 
use the more traditional approach of exactly calculating 
the ground-state properties and the Hamiltonian evolu- 
tion of hardcore bosons for a single realization of disorder 
at a time, averaging then over the disorder distribution 
through Monte Carlo importance sampling. Accepting 
the overhead of disorder averaging, this approach has the 
advantage that, unlike the method of Ref. [H, the time 
evolution is exact for arbitrarily long times. 

According to Eq. ([5]), the weights of the disorder confi- 
gurations are defined by the coefficients of the Fock-state 
decomposition for the initial state of the frozen bosons 



through P({r4}) 
Fock state 



l c ({ n i})l ■ Introducing the bosonic 



and the iV f x L matrix 
some algebra we get 



{<}) = 4A-~ a tJ°) ( 17 ) 



c(K}) = (H}|$ f )=det (Qtp) 



(18) 



The energy eigenstates of noninteracting fermions with 
periodic boundary conditions, contained in the columns 
of the P matrix, are given by plain waves. Thus it can 
be shown that the Slater determinant of Eq. (fTS| is a 
Vandermonde determinant, which can be evaluated ana- 
lytically. The disorder weights finally become 



(HI) I 
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JW 



n sin2 

l<n<m<7V f 



(in 



(19) 



Throughout the paper we will compare the effect of 
the disorder potential generated by the frozen superfluid 
state of N bosons with that of the potential generated 
by N l fully uncorrelated frozen bosons^, with the flat 
distribution 



:(K})| 2 = (iV f !(L-7V f )!)/L! 



(20) 



Even when equipped with the exact statistics of dis- 
order as in Eqs. (fT9]) and (|20|) . it is generally hopeless 
to fully average the results of any calculation on the 
L!/(iV f !(L — N { )\) disorder configurations. Hence we 
opt for Monte Carlo importance sampling, where starting 
from an initially chosen random configuration we 

propose a new one | {^f} ) (e.g. by swapping the occupa- 
tion of two sites) and accept it with Metropolis probabil- 
ity p = min (jc({n-}')| 2 /|c ({"<*}) | 2 , l) • The real-space 

properties are typically averaged over 10 disorder real- 
izations, obtained one from the other by updating O(L) 
sites, while the OPDM is averaged over 10 2 realizations 
(due to the computational overhead of the OPDM cal- 
culation) . This provides full convergence of the disorder- 
averaged quantities. 



D. Characteristics of Disorder 

A fundamental figure of merit for the proposed setup 
to introduce disorder in optical lattices is represented by 
the correlation properties of the disorder potential gen- 
erated by the frozen particles. As discussed in the intro- 
duction, the aspect of correlations represents one of the 
main weaknesses of the optical implementation of disor- 
der through laser speckles 12 . In particular, we focus on 
the density-density correlation function 



C r (|*'» 



f ) (» 



(21) 



where n f = (nl) — N* /L. This correlation function for 
the frozen hardcore bosons is identical to the one for the 
corresponding fermions, and can be calculated exactly in 
the homogeneous system {V 1 — 0) with periodic bound- 
ary conditions; for N f < L, namely in the superfluid 
state of hardcore bosons, 



C r ({&)) = 



for r = 

ItfiL-N*)]- 1 ^T^) forr>l 

. (22) 

Hence, in the limit r <C L (which is always satisfied in the 
thermodynamic limit) the correlator decays algebraically 
like C r (1^-')) ~ r~ 2 with a superimposed oscillation at 
twice the Fermi wavevector kp — 7rn f . The negative val- 
ues of the correlator can be understood as resulting from 
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an effective long-range repulsion between the hardcore 
bosons due to the kinetic energy, which enjoys as much 
free space around each boson as possible. 

From Eq. (|22|) we deduce that the disorder potential 
created by the frozen bosons has slowly decaying cor- 
relations. Nonetheless it has fast oscillations on short 
distances, which are captured by the structure factor, 
namely the Fourier transform of the correlator: 



L-l 



Cfc = ^ C r cos (kr) 



r=0 



This function can be evaluated to give for k £ 
and n l < 0.5 



Cfc 



i(l-n f ) 

l 

L(l-n f ) 



JM. 

2fe F 



for |jfe| < 2k F 

for 2k F < \k\ < it. 



(23) 



-7T, 7r] 



(24) 



The expression for n > 0.5 follows from particle-hole 
symmetry. The corresponding functions for an uncorre- 
cted random potential generated by randomly displayed 
frozen particles are C r = 5 t q and Ck = l/L. We observe 
that the Fourier spectrum of the frozen-boson potential is 
extremely broad, and it has a flat maximum for k > fcp, 
reflecting the short wavelength oscillations of the poten- 
tial. Hence we are in presence of a correlated but strongly 
fluctuating potential. In the following section we will see 
that such a potential generally leads to Anderson local- 
ization of all single-particle states and to the suppression 
of quasi-condensation at all fillings. Throughout the rest 
of the paper, we will take the case of half-filling for the 
frozen bosons, iV f = L/2. 



III. GROUND-STATE PROPERTIES 

In this section we investigate the properties of the 
eigenstates of the Hamiltonian Ho + Hint (see equations 
@ and ([3])) for the mobile bosons moving in the poten- 
tial created by the frozen bosons. Here we consider an 
homogeneous system (V — V { — 0) with periodic bound- 
ary conditions for both species of bosons. We start with 
an analysis of the localization properties in real space 
and continue with an analysis of coherence in momen- 
tum space. 

We investigate the localization properties in the frozen- 
boson potential through the participation ratio (PR) de- 
fined as 



PR (|#)) 



(Eti< 



(25) 



where (rij) is the average particle density on site i. In 
the case of a rectangular density profile the participation 
ratio gives the support of the density profile, for expo- 
nentially localized states it is proportional to the decay 
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FIG. 1: (Color online) Average participation ratio (upper 
plot) and density of states (lower plot) of the single-particle 
energy eigenstates in the correlated random potential created 
by frozen bosons with disorder strength W = 0.5J. The den- 
sity of states is calculated for a system size L = 2002. 



length of the wavefunction, and for a Gaussian-shaped 
density profile it is proportional to its standard devia- 
tion. 

The average PR of single particle eigenstates in 
the potential created by a half-filled system of frozen 
bosons at W = 0.5J is shown in the upper panel 
of Fig. U as a function of energy, for various system 
sizes. The PR essentially becomes size independent 
for L > 1000, clearly indicating the localization of all 
single-particle eigenstates, except for possibly a subset 
of zero measure; we find a similar result for all strengths 
of the potential we investigated. Still a non-extensive 
number of extended eigenstates, not captured by this 
analysis, is in principle sufficient to suppress localization 
e.g. in transport experiments^, so that further analysis 
is required to complete the picture on the localization 
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properties of the frozen-boson potential (see Sec. fV|) . 
The peak in the PR at the center of the band can be 
understood by the fact that the corresponding states 
have dominant k — ±7r/2 components, which are also the 
dominant Fourier components of the random potential 
for half filling. 

We now move on to discuss the properties of the many- 
body ground state for hardcore bosons in the frozen- 
boson potential. Making contact with the discussion of 
Sec. Ill Al such a ground state can be reached by adia- 
batically turning on the interaction W between the two 
species after having prepared each of them separately 
and having quenched the hopping of the frozen bosons. 
In particular we analyze the properties of the disorder- 
averaged OPDM, which is translationally invariant, so 
that its eigenvalues correspond to the momentum distri- 
bution function (MDF), and hence average condensation 
properties are studied in momentum space through the 
scaling of Ao = rifc=o with the number of particles N. 

The scaling analysis of the disorder-averaged occupa- 
tions of the lowest natural orbitals Ao is shown in Fig. [3] 
where, as N grows, N { and L are grown correspond- 
ingly such that N — N l = L/2, For small system sizes 
the hardcore bosons show quasi-condensation behaviour 
Ao ~ N a , with a = 0.5 within the error given by the sim- 
ulation, but for larger system sizes Ao saturates, hence re- 
vealing the absence of quasi-condensation in the thermo- 
dynamic limit. As shown in Fig. [2J the same quantity for 
the potential generated by fully uncorrelated frozen parti- 
cles gives a completely analogous picture. This crossover 
can be qualitatively explained as a fragmentation effect. 
The lowest natural orbital in each disorder configuration 
is localized, namely it does not scale with the system size 
and hence it can can only accommodate a finite number 
of particles, since the bosons repel each other. 

As translational invariance is restored after disorder- 
averaging, the eigenvalues of the disorder-averaged 
OPDM correspond to the momentum distribution func- 
tion. The peak at zero momentum, in which quasi- 
condensation appears in the case without disorder, is 
significantly reduced in presence of disorder. A scaling 
analysis of nk=o of the mobile bosons, shown in Fig.[5]de- 
tects a similar crossover from algebraic increase to satura- 
tion as for the disorder-averaged Ao- Therefore, the con- 
sequences of fragmentation extends to this experimen- 
tally accessible observable, and we can conclude on the 
absence of quasi-condensation in the disorder-averaged 
OPDM of the ground state of the system^. Moreover 
the density of states, shown in the lower panel of Fig. [1] 
reveals a continuous excitation spectrum at all energies. 
Hence, for a system of interacting bosons, the absence 
of quasi-condensation, together with the absence of gaps 
in the density of states, leads to classify the state of the 
system as a Bose glass. At the same time the Jordan- 
Wigner transformation translates this phase into an ideal 
Anderson insulator of non- interacting fermions. 




1 J , r 

10 100 
N 

FIG. 2: (Color online) Scaling of the disorder-averaged largest 
eigenvalue Ao of the OPDM (upper plot) and the occupation 
of the zero momentum state rik=o (lower plot) for the ground 
state of a half-filled system of mobile bosons, interacting with 
the disorder potential. The cases of a correlated potential 
resulting from the superfluid state, Eq. (|19|l . and that of a 
fully uncorrelated potential, Eq. (|20|l . are compared. Here 
the strength of disorder is W = 0.5 J. The black lines are fits 
to rik=a, (Ao) oc i/N for the first four data points. 



IV. DYNAMICAL PROPERTIES AFTER 
SUDDEN ON-TURN OF THE DISORDER 
POTENTIAL 

While the previous section investigated the state of the 
mobile bosons after an adiabatic on-turn of the interac- 
tion with the disorder potential created by the frozen 
bosons, in this section we consider the evolution of the 
mobile bosons after turning on that interaction suddenly. 
We consider such a time evolution in the homogeneous 
system, i.e. V = V 1 = 0, with periodic boundary con- 
ditions, mimicking the situation in which both species 
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FIG. 3: (Color online) Scaling of the steady-state zero- 
momentum occupation nt=o with the particle number TV in a 
homogeneous system of mobile and frozen bosons (V = V f = 
0) at half-filling (N = iV^ = L/2), with interaction strength 
W — 0.5J. The line corresponds to a fit nk=o ~ y/N to the 
first four data points. 



FIG. 4: (Color online) Snapshot of the evolution of the real- 
space density rn for a single particle, initially in the ground 
state of a harmonic trap with V = 0.01J, and interacting with 
a correlated disorder potential. Here the interaction strength 
is W — 0.5J, and the system size is L = 502. 



of particles are prepared in the same region of space be- 
fore being brought into interaction. In this situation, the 
disorder-averaged real-space density of the mobile bosons 
remains uniform during time evolution, so that real-space 
localization effects are not visible. Nonetheless, even in 
a more realistic experimental scenario, in which both 
species of particles are kept in the same nonvanishing 
trap (V — V i > 0), the disorder-averaged density profile 
of the mobile species does not reveal marked localization 
effects, and it is even expanding during time evolution 
in order to reduce the overlap with the confined frozen 
boson a 28 ! 44 . 

The fundamental effect of disorder on the time evolu- 
tion can be very clearly detected in momentum space, 
which can be measured in time-of- flight experiments. 
The initial quasi-condensation peak in the MDF at k — 
decreases quickly during time evolution until a stationary 
regime is reached, in which the value of rik=o is oscillating 
with a small amplitude (~ 1%). The scaling analysis for 
the value of nk=o, time-averaged over the small oscilla- 
tions, is shown in Fig. [3] for a half-filled system of mobile 
bosons. The occupation of the zero momentum state 
n/c=o increases approximately like y/~N for small N, but 
it then deviates from an algebraic increase and saturates 
for larger N; a completely analogous behavior is found 
for the case of uncorrelated disorder. Hence the scaling 
of rik=o reveals a crossover from quasi-condensation to 
fragmentation, with close similarity to the case of adia- 
batic evolution described in Sec. IIIII in analogy to that 
case, we conclude that the steady state reached by time 
evolution realizes dynamically a Bose glass. 



V. TRANSPORT PROPERTIES 

In this section we discuss the expansion properties of 
initially confined mobile bosons in the potential created 
by the frozen particles. The initial confinement for the 
mobile particles is provided by a tight parabolic trap 
(V > 0), while the frozen particles can be imagined as 
prepared in a much shallower trap, whose effect is ignored 
for simplicity, so that the random potential they gener- 
ate is the same as the one studied in the homogeneous 
case in Sec. Ill Dl In particular we imagine that at time 
t = the trap confining the mobile bosons is released 
(V = 0) and simultaneously the two species are brought 
into interaction (W > 0). After release from the trap, 
the mobile particles would expand forever in absence of 
disorder, so that a halt in presence of disorder explicitly 
shows localization effects. 

Subsection IV Al deals with the evolution of the real- 
space properties during expansion, while Subsection IV Bl 
studies the evolution of the condensation properties of 
the system. 



A. Real-space localization 

In this subsection we focus on the real-space properties 
of the time-evolved hardcore bosons in the frozen-boson 
potential, starting from a confined state in a parabolic 
trap. Snapshots from a disorder-averaged time evolution 
of a single particle starting from the ground state in the 
trap are depicted in Fig.2) It is shown that the expansion 
reaches a localized steady state in the long-time limit; in 
particular the spatial decay of the particle density under- 
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FIG. 5: (Color online) Time evolution of the participation 
ratio of a single particle for various disorder strengths. Other 
parameters as in Fig. [4] 



goes a dramatic change from Gaussian - as expected in a 
parabolic trap - to exponential - as expected in presence 
of Anderson localization. Hence an Anderson localized 
state is realized dynamically during expansion; it is fun- 
damental to stress that this steady state does not corre- 
spond to the ground state of the system, as the energy 
of the particle is conserved during expansion and hence 
it does not relax to the ground-state value. 

The dynamical localization of the single particle wave- 
function is fully captured by the time evolution of the PR, 
depicted in Fig. [5] for various strengths W of the disorder 
potential. Without disorder, the wavefunction spreads 
ballistically without changing its Gaussian shaped. In 
the presence of disorder, the PR saturates instead to a 
finite value, which decreases when increasing the disor- 
der potential; a fit to the final value of the PR, shown in 
Fig. [Bl suggests that saturation in the PR takes place for 
any arbitrarily small value of the potential, as it would 
be expected for uncorrelated disorder, although explor- 
ing very small strengths of the disorder is numerically de- 
manding as the steady state is reached for prohibitively 
large values of the PR. 

The time evolution of the participation ratio of many 
particles initially confined in a harmonic trap is similar 
to that exhibited by a single particle, and a saturation 
of the PR to a steady-state value is observed for the 
smallest value of W that we could treat numerically 
(Fig.©. 

Nonetheless, upon increasing the particle number in 
the trap beyond a critical value N c « 26 (see equation 
[14] with V = 0.01) a Mott plateau at unit filling starts 
forming in the center, and the properties of the initial 
trapped ground state change drastically. It is then inter- 
esting to study whether the different initial conditions for 



FIG. 6: (Color online) Steady state participation ratio 
(recorded at t = 500J - ) for N particles, time-evolved in 
a system of size L = 702 starting from the ground state of 
a harmonic trap (V = 0.01J) and interacting with a corre- 
lated random potential or with a fully uncorrelated random 
potential. The lines show a fit to the data points with an 
algebraically decaying function. 



the time evolution reflect themselves in the steady state 
of the system after the halt of the expansion. The upper 
plot of Fig. [7] shows the PR of the steady state as a func- 
tion of the number of particles N. The steady state PR 
increases with the number of particles for N < N c , show- 
ing a small peak for N — N c beyond which the increase 
with N is much slower. This feature can be understood 
in terms of the time evolution of the fermionic wavefunc- 
tion Eq. (fTTj) , whose real-space properties such as the PR 
are equivalent to those of the hardcore bosons, and whose 
time evolution is simply obtained through the time evo- 
lution of single-particle eigenfunctions. For N > N c the 
single particle wavefunctions that are populated for in- 
creasing N are more and more confined to the sides of the 
trap, as the center has already a saturated density, and 
their energy is dominated by the trapping term, while the 
single-particle kinetic energy decreases with increasing N 
beyond N c (see lower plot of Fig. [7]). Given that the trap- 
ping energy vanishes at t = after the trap release, the 
fermions occupying levels beyond the iV c -th expand with 
an energy which is less than that of the iV c -th level, so 
that they are expected to be localized by the disorder 
potential with a final participation ratio which is similar 
(or even less than) that of the iV c -th particle. The time- 
evolved particle density of the many-body system is the 
simple sum of the squared time-evolved wavefunctions of 
the free fermions, so that the total participation ratio is 
not expected to increase significantly when adding parti- 
cles beyond the iV c -th. 
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FIG. 7: (Color online) Upper plot: Steady-state participation 
ratio (recorded at t = 400 J -1 ) of N particles as a function 
of the particle number for various disorder strengths. Other 
parameters as in Fig. [6] Lower plot: Total, kinetic, and poten- 
tial energy of the Mth single-particle eigenstate of the initial 
system in the trap. 



B. Coherence properties 

In the previous section the real-space properties of ini- 
tially confined hardcore bosons expanding in a disorder 
potential have been discussed. This section addresses 
the condensation properties of the same system, moti- 
vated by the rich physical scenario offered by expanding 
hardcore bosons in absence of disorder— 122. Subsection 
IV B H is devoted to the time evolution of hardcore bosons 
initially confined in a superfluid state, while Subsection 
IV B 21 analyses the expansion from a Mott insulator. 



1. Time evolution starting from a superfluid state 

This subsection focuses on the expansion of the hard- 
core bosons initially confined in a shallow trap with char- 



FIG. 8: (Color online) Snapshots of the time evolution of the 
momentum distribution rik for a system oi N — 100 hardcore 
bosons, initially in a trap with density p — 0.5, interacting 
with a correlated disorder potential. Here the system size is 
L = 1002 and the interaction strength is W = 0.5J. The mo- 
mentum distribution for the corresponding spinless fermions 
is shown for comparison. 



acteristic density p = 0.5 (see Eq. (|14l) ) well below the 
critical value p c » 2.6 for the onset of a Mott plateau 
in the trap center. With this initial condition quasi- 
condensation in the first NO is present, and the MDF 
is peaked around k = at t = 37 . In the absence of 
disorder the MDF evolves towards that of the fermions 
(which is a constant of motion), namely the bosons 
fermionize in momentum space, although they still quasi- 
condense in the lowest NO which contains many momen- 
tum components 3 ^. Furthermore, the OPDM is decaying 
algebraically like pi^+ r ~ M -0 ' 5 , showing quasi-long- 
range order, whereas the phase of the OPDM is oscil- 
lating at large distances leading to fermionization in mo- 
mentum space. The decay of the OPDM is measured 
from the center of the trap, i.e. the center of the lowest 
NO and the site with maximum occupation. 

We now move on to the analysis of the evolution of 
the MDF during expansion from a shallow trap in the 
presence of disorder created by the frozen bosons. The 
MDF of the mobile hardcore bosons is depicted in Fig. [8] 
at different times t, compared to the MDF of fermions 
in the same system. In presence of disorder the MDF 
of the fermions is no longer a constant of motion, but 
it is only slightly broadening in time due to its interac- 
tion with the frozen species of particles. The effect of 
the disorder potential on the hardcore bosons is more 
significant. The t = peak at zero momentum reduces 
its height, but, contrary to the expansion without dis- 
order, it does not disappear and fermionization is not 
present. Following Ref. 1371 , fermionization in absence of 
disorder is understood through the argument that, after 
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FIG. 9: (Color online) Decay of the modulus of the OPDM 
|/Oj ,i +r| from the system center io at various times during 
evolution. Parameters as in Fig. [8] 



FIG. 10: (Color online) Scaling of the steady-state Ao 
(recorded at time t = 500 J ) for varying number of par- 
ticles in a system of size L = 1502; other parameters as in 
Fig. [8] The data for correlated frozen bosons are compared 
with those for fully uncorrelated ones. The line corresponds 
to a fit to Ao ~ vN. 



a long-term expansion, the hardcore-boson system is di- 
lute enough to be considered essentially non-interacting, 
and so it becomes equivalent to its fermionic counterpart. 
In presence of disorder, on the contrary, the expansion 
is stopped by localization, so that the extremely dilute 
limit is never reached and the hardcore bosons preserve 
their nature of strongly interacting particles. Further- 
more, the interaction with a disorder potential leads to 
the loss of quasi-long-range order during time evolution, 
as shown in Fig. [9] by the decay of the OPDM at various 
times. The system finally reaches a steady state with an 
exponentially decaying OPDM. 

The loss of quasi-long-range order during expansion 
strongly suggests the loss of quasi-condensation in the 
system. As the initial conditions break the translational 
symmetry, condensation properties are not captured by 
the scaling of the occupation of the zero-momentum 
state, and direct diagonalization of the disorder-averaged 
OPDM is necessary to extract the scaling of the occupa- 
tion Ao of the lowest natural orbital. Ao decreases during 
time evolution, reaching a constant value in correspon- 
dence with the steady state observed in real-space. The 
scaling analysis of Ao is performed for this constant value 
in Fig. [TO] The scaling of Ao deviates from the quasi- 
condensation behavior at sufficiently large particle num- 
bers, revealing a crossover from quasi-condensation to 
fragmentation due to the localized nature of the lowest 
NO. Fully uncorrelated disorder leads to a similar behav- 
ior, although the crossover appears to be much broader, 
and a quasi-condensation regime at low N could not be 
identified. 



2. Time Evolution starting from a Mott insulator 

As reported in Ref. HB, hardcore bosons prepared in 
a perfect Mott insulating state (corresponding to an in- 
finitely steep trap) and subsequently time evolved exhibit 
the phenomenon of dynamical quasi-condensation at fi- 
nite momentum. Indeed, from the initially flat MDF of 
the Mott insulator two quasi-condensation peaks emerge 
at momenta k = ±7r/2 during time evolution— , as repro- 
duced in Fig. [TTJ Strictly speaking quasi-condensation 
happens in two degenerate NOs whose Fourier transform 
is sharply peaked around k = ±7r/2; the NOs propagate 
at a velocity v = ±2 J corresponding to the maximal 
group velocities dek/dk for the single-particle dispersion 
relation Ck — — 2 J cos k at momenta k = ±7r/2. 

The occupation of the degenerate lowest NOs, Ao, fol- 
lows initially a universal power-law increase with time, 
independent of the number of particles N; at an N- 
dependent characteristic time r c , the two degenerate low- 
est NOs begin to move in opposite directions, and Ao 
starts to algebraically decrease, although its scaling with 
particle number shows the typical quasi-condensation be- 
havior Ao ~ y/N. The two lowest NOs clearly appear in 
the real-space densities as coherent fronts of the atomic 
cloud moving in opposite directions (see Fig. ITT]) . This 
aspect suggests the use of this setup to produce an atom 
laser—. 

We now consider the case of expanding hardcore 
bosons in the disorder potential created by frozen par- 
ticles. Fig. [T2] shows the density profile and the MDF 
for the expansion from a perfect Mott insulator in the 
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FIG. 11: (Color online) Time evolution of the density pro- 
file (upper image) and the momentum distribution function 
(lower image) of N — 35 hardcore bosons, initially in a perfect 
Mott insulator, in a system of size L — 502 without disorder. 



FIG. 12: (Color online) Time evolution of the density pro- 
file (upper image) and the momentum distribution function 
(lower image) of N = 35 hardcore bosons, initially in a Mott 
insulator, interacting with the disorder potential created by 
frozen particles at half-filling. Here the interaction strength 
is W = 0.5J. 



presence of disorder created by frozen particles. Initially 
peaks at momenta k = ±7r/2 are emerging from the flat 
MDF at t = as in the case without disorder, and the 
hardcore boson cloud shows two outer fronts expanding 
ballistically in opposite directions. Nonetheless this ini- 
tially coherent expansion is rapidly suppressed due to 
localization, which leads to a decrease and broadening of 
the momentum peaks at larger times up to a final steady- 
state MDF in which two broad peaks survive, and which 
corresponds to a localized state in real space as seen in 
Subsection IV Al 

A deeper analysis of condensation effects relies on the 
NOs and their occupations. Figs. [T3] and [T3] show the 
time evolution of the largest eigenvalue Ao for different 
disorder strengths and different particle numbers. Ini- 
tially Ao is two-fold degenerate, corresponding to reflec- 
tion symmetry at the center of the system, and it in- 
creases following a universal power law independent of 
the particle number, similarly to what is observed in ab- 
sence of disorder but with a different disorder-dependent 
exponent. As in the case W — 0, at a characteristic time 
t c the two degenerate natural orbitals start to move into 
opposite directions, and correspondingly the time evolu- 
tion of Ao turns into a decreasing behavior; r c strongly 
depends on the disorder strength (Fig. I13p. while its de- 
pendence on the particle number N becomes weaker for 



large N, where r c is seen to approach an asymptotic value 
(Fig. rnT| . Unlike the case W = 0, at a later stage the 
expansion of the system is stopped by disorder, and the 
degeneracy in the lowest NOs is removed, going from two 
propagating ones to a single NO localized in the system 
center. 

Hence the expanding system of hardcore bosons from 
a pure Mott state shows a crossover from an incipient 
quasi-condensation regime at finite momenta to a local- 
ization regime. It is natural to ask whether the sys- 
tem displays true quasi-condensation at any intermediate 
point in time. To address this issue, we perform a scal- 
ing analysis of Ao at the maximum-coherence time t = t c . 
The results of this analysis are shown in Fig. 1151 where 
Ao exhibits a clear saturation for large particle numbers, 
and hence a crossover from quasi-condensation to frag- 
mentation. Repeating the same analysis for uncorrelated 
random disorder we find a scaling of Ao which is in qual- 
itative agreement with the case of disorder generated by 
frozen particles (the apparent quantitative agreement is 
accidental) . 
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FIG. 13: (Color online) Time evolution of Ao for various dis- 
order strengths. Other parameters as in Fig. 1121 



FIG. 15: (Color online) Scaling of the maximal value of Ao 
during time evolution of N hardcore bosons, initially in a 
Mott insulator. All parameters as in Fig. [14] . The line corre- 
sponds to a fit to Ao ~ ^iV of the first four data points. 
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FIG. 14: (Color online) Time evolution of Ao for various num- 
bers of particles, and for disorder strength W = 0.1J. Other 
parameters as in Fig. 1121 



VI. EXPERIMENTAL REALIZATION 

The original experimental proposal motivating this 
work is reported in Ref. [H, where the use of state- 
dependent lattices was initially envisioned to create the 
quantum superposition of random potentials. As shown 
in Ref. 31, two hyperfine states of the same atomic species 
can be loaded in the minima of two different polarization 
components of an optical lattice at a so called "magic 
wavelength" 5 i 29 i 30 , at which each hyperfine state couples 
to one and only one polarization component. Hence, if 
the two polarization components are initially shifted by 
it/2 in space, the two species are non essentially inter- 
acting (see Fig. \W\i . The possibility of increasing dras- 



tically the intensity of one of the two polarization com- 
ponents would allow for a sudden quench of the hop- 
ping of one of the two species, preparing in this way the 
quantum superposition of random potentials. The two 
species can be then brought into interaction at differ- 
ent strengths by shifting the spatial phase between the 
two polarization components of the lattice. An adiabatic 
shift would transfer the mobile species to the Bose-glass 
ground state in the random potential, as discussed in 
Section IIII1 whereas a sudden shift would give rise to 
an out-of-equilibrium Bose-glass state as steady state af- 
ter a transient evolution, as discussed in Section [TVl As 
demonstrated in Ref. l33l . the hardcore regime is achieved 
by using deep lattices for both polarization components 
and extremely dilute gases of both hyperfine states, so 
that a density of less than one atom per site is achieved 
in the lattice after loading. 

As seen in Section [Vj expansion experiments in the 
disorder potential require to confine the two components 
with strongly different trapping frequencies, and then to 
release one of the two traps independently of the other. 
Making use of the selective coupling of different hyperfine 
states to different polarization components of a magic- 
wavelength laser, it is possible to more strongly con- 
fine one of the two species through an optical dipole 
trap obtained by a tightly focused and polarized laser 
propagating transverse to the chain direction in the one- 
dimensional optical lattice. 

Throughout all the previous Sections we have seen 
that the onset of localization leaves very strong signa- 
tures on the momentum distribution, which is observed in 
standard time-of- flight experiments. Moreover expansion 
experiments in the disorder potential lead to exponen- 
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FIG. 16: (Color online) Sketch of the preparation of the dis- 
order potential and of the tuning of its strength using state- 
dependent optical lattices for a mixture of atoms in two inter- 
nal states. The two species are first prepared in two shifted 
optical lattices (a), then one of the two optical lattice is in- 
creased in strength so as to freeze one of the two bosonic 
species (b), and finally the spatial phase of the two lattices is 
changed to bring the two species into interaction (c). 

tially localized steady states, whose density profile can 
be reconstructed by absorption images taken shortly af- 
ter turning off all lasersl&i^. Despite the averaging over 
the various tubes of the one-dimensional optical lattice 33 , 
the extreme tails of the averaged particle density distri- 
bution should be dominated by the exponentially local- 
ized tails of the particle density in the central tubes of 
the one-dimensional optical lattice. 

VII. CONCLUSIONS 

In this paper we have shown that a gas of one- 
dimensional hardcore bosons undergoes genuine quantum 
localization effects when set into interaction with a sec- 
ondary species of bosons frozen in a massive quantum 
superposition of Fock states. Each Fock state can be re- 
garded as a realization of a random potential, and the 
unitary evolution of the mobile species of bosons follows 
all possible paths related to the various disorder realiza- 
tions in parallel. Physically relevant states in which one 



can prepare the frozen bosons, as for instance the super- 
fluid state in the hardcore limit, realize a rapidly fluctuat- 
ing disorder potential over the length scale of a few lattice 
spacings; despite its power-law decaying correlations, the 
disorder potential is found to lead to the same localiza- 
tion effects as those observed in a fully uncorrelated po- 
tential. In the hardcore boson limit for the mobile species 
these effects can be studied exactly in real time through 
Jordan- Wigner diagonalization, and we can numerically 
simulate realistic experimental setups for the dynamical 
preparation of localized many-body states. The equilib- 
rium state of the hardcore bosons in the random poten- 
tial is found to be a homogeneous Bose-glass state with 
exponentially decaying correlations; a similar state can 
be realized also dynamically after a sudden on-turn of 
the interaction between the two species. When the hard- 
core bosons are initially confined in a tight trap and then 
set free to expand in the random potential, for any non- 
vanishing disorder strengths the expansion stops and the 
system reaches an exponentially localized state. 

For any setup discussed in this work, we observe the 
absence of quasi-condensation and quasi-long-range or- 
der due to the disorder potential. When the two species 
are confined in the same region and brought into inter- 
action, the steady state of the system displays exponen- 
tially decaying off-diagonal correlations; in the expansion 
experiments the disorder potential destroys the effects of 
fermionization (when expanding from an initially dilute 
state) and quasi-condensation at finite momentum (when 
expanding from a Mott state). Hence disorder created 
by a species of frozen hardcore bosons represents a very 
robust way to experimentally implement strongly fluctu- 
ating random potentials in optical lattices, and to realize 
fundamental localization effects of many-body systems. 
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